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c3 1. Introduction 



ABSTRACT 



We performed a uniform and detailed analysis of 1 1 12 high-resolution spectra of FGK dwarfs obtained with the HARPS spectrograph 
at the ESO 3.6 m telescope (La Silla, Chile). Most stars have effective temperatures 4700 K < T eff < 6300 K and lie in the metallicity 
range of -1.39 < [Fe/H] < 0.55. Our main goal is to investigate whether there are any differences between the elemental abundance 
trends (especially [a7Fe] ratio) for stars of different subpopulations. The equivalent widths of spectral lines are automatically measured 
from HARPS spectra with the ARES code. The abundances of three a elements are determined using a differential LTE analysis 
relative to the Sun, with the 2010 revised version of the spectral synthesis code MOOG and a grid of Kurucz ATLAS9 atmospheres. The 
stars of our sample fall into two populations, clearly separated in terms of [a/Fe] up to super-solar metallicities. In turn, high-or stars 
are also separated into two families with a gap in both [a/Fe] ( [a/Fe] « 0.17) and metallicity ([Fe/H] = -0.2) distributions. The 
metal-poor high-a stars (thick disk) and metal-rich high-a stars are on average older than chemically defined thin disk stars (low-a 
stars). The two a-enhanced families have different kinematics and orbital parameters. The metal-rich ff-enhanced stars, such as thin 
disk stars have nearly circular orbits, close to the Galactic plane. We put forward the idea that these stars may have been formed in 
the inner Galactic disk, but their exact nature still remains to be clarified. 

Key words, stars: abundances - stars: kinematics and dynamics - Galaxy: disk 
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The investigation of stellar populations is very important to un- 
derstand the formation and evolution of our Galaxy. The Milky 
Way (MW) has a composite structure with several subsystems. 
The main three stellar populations of the MW in the solar neigh- 
borhood are the thin disk, thick disk, and the halo. These pop- 
ulations have different kinematic and chemical properties. The 
subdivision between the disk and halo was first identified since 
long ago, but the thick disk was discovered far more recently by 
Gilmore & Reid (1983 ), who analysed the stellar density distri- 
bution as a function of distance from the Galactic plane. 

There is no obvious predetermined way to identify purely 
thick or thin disk stars in the solar neighborhood. There are 
essentially three ways of distinguishing local thick and thin 
disk stars: a purely kinematical approach (e.g. Bensby et al. 
120031 [20531 Reddy et a l. |2006), a purely chemical method 
(e.g. Navarro et al. 1201 II) . and by looking at a combination of 
kinematics, metallicities, and stellar ages (e.g. Fuhrmann 1998; 
Haywood|2008il). 

The kinematic selection is a much more commonly applied 
method than the chemical approach, because it is much easier 
to measure the velocity of a star than to determine its chemi- 
cal composition (particularly its ^-enhancement). However, the 
chemical distinction of the disks can be more useful and reliable, 
at least, because chemistry is a relatively more stable property 
of a star, that is intimately connected to the time and place of 
its birth, whereas spatial positions and kinematics are evolving 
properties. 



During the past few years, there have been several studies 
directed to the detailed elemental abundance investigations of 
stars in different subpopulations. However, spectroscopic stud- 
ies are in general limited to small samples of a few hundred 
stars at most (e.g. Feltzing & Gustafsson 1998; Bensby et al. 
120031 120051 120071 Soubiran et al. 120051 Reddy et al. 120061 
Ramirez et al. 2007) and only a few studies have been based 
on samples as large as 1000 stars (e.g. Gazzano et al. 2010; 
Gazzano 120111 Petigura & Marcy |201 II ). To investigate [a/Fe] 
abundances in the thin and thick disks with relatively large sam- 
ples, some studies combine data from different sources (e.g. 
Navarro et al. 1201 II ). Alternatively, some other studies estimated 
the [a/Fe] ratio from Stromgren indices (e.g. Casagrande et al. 
1201 II) . However, both methods are far less precise than those 
obtained with high-resolution spectroscopy, and prevent us from 
seeing any separation gap between the thin and thick disks. To 
minimize any type of external and internal "errors", one needs 
to have as large and as homogeneous a sample as possible, with 
reliable measurements of their chemical and kinematic features. 

In this Letter, we investigate the possible differences in the 
elemental abundance trends for stars of different subpopulations, 
using a stellar sample of 1112 long-lived dwarf stars. To sepa- 
rate and investigate the different Galactic stellar subsystems, we 
focus on the [a/Fe] ratio (here "a" refers to the average abun- 
dance of Mg, Si, and Ti). The extensive and full investigation of 
this sample, will be more focused on the abundance difference 
between stars with and without planets and be presented in an 
upcoming paper where we will also describe the observations, 
data reductions, and abundance analysis in detail. 
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2. The sample 

The sample used in this work consists of 1112 FGK stars ob- 
served within the context of the HARPS GTO programs, here- 
after called HARPS- 1 (Mayor et al. l2003l ), HARP S-2 (Lo Curto 
et al. l20T0l l. and HARPS-4 (Santos et al. [2QTTb . The stars are 
slowly-rotating, non-evolved, and in general have a low level of 
activity. The individual spectra of each star were reduced using 
the HARPS pipeline and then combined with IRAF0 after cor- 
recting for its radial velocity. The final spectra have a resolution 
of R ~l 10000 and signal-to-noise ratio (S/N) ranging from ~20 
to ~2000, depending on the amount and quality of the original 
spectra. Fifty-five percent of the spectra have S/N higher than 
200. 

Precise stellar parameters for all the stars were determined 
in the same manner and by the same authors from the same 
spectra used in our study. For details, we refer to Sousa et al. 
(2008] 1201 lal 1201 lb\ . The typical uncertainties in the atmo- 
spheric parameters are of the order of 30 K for T^g , 0.06 dex 
for log g, and 0.03 dex for [Fe/H]. 

The stars in the sample have effective temperatures 4487 K 

< T eS < 7212 K (there are only 12 stars with T eS > 6500 K) and 
metallicites -1.39 < [Fe/H] < 0.55 (only 11 stars with [Fe/H] 

< -1 and three stars with [Fe/H] > 0.4), and they have surface 
gravities 2.68 < log g < 4.96 dex (again the number of "outliers" 
is very small, only 5 stars with log g < 3.8 dex). 

3. a element abundance 

Elemental abundances for three a elements (Mg, Si, and Ti) were 
determined using a differential LTE analysis relative to the Sun 
(the reference abundances were taken from Anders & Grevesse 
dl989l l) with the 2010 revised version of the spectral synthesis 
code MOOGpkSneden [19731 and a grid of Kurucz ATLAS9 at- 
mospheres (Kurucz et al. 1993). The method and the atomic data 
are the same as in Neves et al. ( 120091 ) (see also Adibekyan et al. 
1201 II in preparation). The equivalent widths (EW) were auto- 
matically measured with the ARES0 code (Automatic Routine 
for line Equivalent widths in stellar Spectra - Sousa et al. 120071 ) . 
The input parameters for ARES were calculated following the 
procedure discussed in Sousa et al. d201 lal) . 

The total uncertainties in the abundances were evaluated as 
combinations of the line-to-line scatter errors and errors induced 
by uncertainties in the model atmosphere parameters. The line- 
to-line scatter errors were estimated as cr/ y/N, where cr is the 
standard deviation of N measurements. The abundances of the 
elements were determined using 16, 3, 24, and 6 lines for Si, Mg, 
Til, and Till, respectively. The average error in the [a/Fe] ratio is 
0.05 dex. We note that, in general, our derived abundances agree 
with the [a/Fe] presented in the literature (see Neves et al. 2009 ; 
Adibekyan et al. 1201 II in preparation). 

Studying the [a/Fe] error distribution with stellar parame- 
ters, we observe no discernible trends except in the case of the 
Jeff where we found that for cooler stars (T e ff < 4900 K) the er- 
ror increases, reaching about 0. 1 dex for stars with T e ff = 4500 
K. We also observe a trend with T e s for [Til/Fe] ratio and we 
decided to establish a cutoff temperature - r cuto ff=4900 K (for 
details see Neves et al. 2009). The number of stars in the sample 
with Jeff > 4900 K is 940. 

1 IRAF is distributed by National Optical Astronomy Observatories, 
operated by the Association of Universities for Research in Astronomy, 
Inc., under contract with the National Science Foundation, U.S.A. 

2 http://www.as.utexas.edu/~chris/moog.html 

3 http://www.astro.up.pt/sousasag/ares 
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Fig. 1. [a/Fe] versus [Fe/H] for the whole sample (middle) and 
for stars with T e g -T Q ±300 K (top). Magenta squares refer 
to halo stars. Blue filled circles are the separation points be- 
tween low- and high-a stars, which are minimas of the [a/Fe] 
histograms for five metallicity bins (from [Fe/H] = -0.7 to 0.25) 
and the blue dashed curve is the corresponding separation curve 
passing on that points. The bottom panel is the [a/Fe] versus 
[Fe/H] plot for the whole sample in several bins of metallicity. 
The blue triangles refer to the thick disk stars, black filled circles 
to the metal-rich high-a stars (hamr), and the red crosses to the 
thin disk stars. Error bars in the lower panel correspond to the 
standard deviations, and the error bar in the upper and middle 
panels are the average errors in the [a/Fe] and [Fe/H]. 
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Fig. 2. High-a and low-a separation histograms for the stars with 
T e s -T & ±300 K after subtracting the separation curve. 

Fig. 1 shows [a/Fe] versus [Fe/H] for all stars in the sam- 
ple and for the stars with effective temperatures close to the Sun 
by ±300 K (the number of these stars is 483). As can be seen, 
the stars are clearly separated into two groups according to the 
content of a elements: the "high-a" and the "low-a" stars (thin 
dislfl). This separation highlights the well-known a enhance- 
ment of thick disk stars relative to the thin disk found for stars 
with [Fe/H] < (e.g. Fuhrmann [T998l Bensby et al. l2003l[2005l . 
The blue filled circles in Fig. 1 are the separation points between 
low- and high-a stars. A separation into these two "populations" 
was performed based on [a/Fe], for the stars with r e ff =T G ±300 
K. We divided the sample into five metallicity bins from [Fe/H] 
= -0.7 to 0.25 (see Fig. 2. for bin sizes) and identified the min- 
ima in the [a/Fe] histograms for each bin. The separation curve 
in Fig. 1 is the simple connection of the above-mentioned sepa- 
ration points. In the metallicity region [Fe/H] < -0.7, all the stars, 

4 We emphasize that the traditional selection (and therefore defini- 
tions) of thin and thick disk stars was based on their kinematics. 
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Table 1. The number of stars in different populations. 
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Fig. 3. The [Fe/H] and [a/Fe] separation histograms for the a- 
enhanced stars. 
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Fig. 4. The Toomre diagram for the entire sample. The blue tri- 
angles and red crosses refer to the chemically selected thick 
and thin disk stars, and the black filled circles are hamr stars. 
Magenta squares refer to halo stars selected by their kinematics. 



except five halo stars, are a-enhanced and lie above the [a/Fe] 
= 0.19 line. A distinction between low-a and high-a stars in the 
supersolar metallicity region ([Fe/H] > 0.25) was made using 
the extrapolation of the separation line constructed from the last 
two separation points. After subtracting the separation curve, we 
again compiled the separation histograms presented in Fig. 2. As 
one can see all the separation minimums are very close to [a/Fe] 
= 0, which proves that the separation is clear. 

It is interesting to see that high-a stars are also divided into 
two subgroups: high-a metal-rich stars (hereafter "hamr"), and 
high-a metal-poor stars (thick disk). Fig. 3 shows the [Fe/H] and 
[a/Fe] separation histograms for the or-enhanced stars. In this 
figure, one can see that there is a gap between these two groups 
in terms of both the metallicity (Fe/H ~ -0.2) and [a/Fe] ([a/Fe] 
~ 0.17 ). This gap (in [a/Fe]) can also be seen in Reddy et al. 
(2006) (Fig. 20), although they differentiated between thin and 
thick disk stars using their kinematics. 

On the basis of the above discussed histograms, we can de- 
fine the metallicity ranges for each population. We find that for 
thick disk stars, metallicities vary from -1.4 to -0.2 dex. For thin 
disk stars, the metallicities range roughly from -0.7 to +0.4 dex. 
The hamr group stars lie in the metallicities range roughly from - 
0.2 to +0.3 dex. These values, in particular the upper limits to the 
thick disk and hamr, and the lower limits to the thin and harm- 
groups, may represent the metallicity boundaries of the different 
populations. 

The magenta squares in Fig. 1 refer to stars belonging to halo 
(selected by their kinematics, see Section 4). As is clearly evi- 
dent, halo stars can also be divided into two high-a and low-a 
groups. This dichotomy confirms the results provided in Nissen 
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& Schuster (2010), where they present clear evidence of two dis- 
tinct halo populations. 

4. Kinematics 

To study the kinematic relations and differences between the 
three above separated groups, we computed their Galactic ve- 
locities. The Galactic space velocity components (UVW) were 
derived with respect to the local standard of rest (LSR), adopt- 
ing the standard solar motion ( U , V Q W ) = (10.00, 5.25, 7.17) 
km s _1 of Dehnen & Binney dl9981 l. The main source of the 
parallaxes and proper motions were the updated version of the 
Hipparcos catalog (van Leeuwen 2007) and the radial velocities 
obtained from the HARPS spectra. The parallaxes with errors 
larger than 10%, (which is true for less than 5% of the stars in the 
sample) were redetermined following the procedure described in 
Sousa et al. (1201 lab . The number of stars with inaccurate proper 
motions (errors larger than 10%) is less than 8%. We did not 
perform a quality selection of them, because these errors in gen- 
eral do not change their memberships of a certain population. 
Combining the measurement errors in the parallaxes, proper mo- 
tions, and radial velocities, the resulting average errors in the U, 
V, and W velocities are 2, 2.3, and 1.5 km s , respectively. 

We adopted the method of assigning the probability of each 
star to either the thin disk, the thick disk, or the halo described 
in Reddy et al. (12006b . This assumes that the sample is a mix- 
ture of the three populations and each population follows a 
Gaussian distribution of random velocities in each component 
(Schwarzschild II 9071) . In this paper, we adopt the mean val- 
ues (asymmetric drift) and dispersion in the Gaussian distribu- 
tion, and the population fractions were taken from Bensby et al. 
(120031) and Robin et al. (120051 ) . We consider that a probability in 
excess of 70% suffices to assign a star to the concrete population. 

These kinematic separation criteria suggest that most of the 
stars in the hamr stellar family, such as chemically defined thin 
disk stars, have thin disk kinematics (see Table 1). The same 
impression is created by Fig. 4, where all the stars in this study 
(with r e ff > 4900 K) are shown in a Toomre diagram. The chem- 
ically defined "thick" disk stars mainly have thick disk kinemat- 



5. Discussion 

The near constancy and the high level of a-element abundances 
relative to Fe for the metal -poor thick disk stars ([Fe/H] < -0.3) 
suggest that they formed in regions of high star formation rate 
where massive stars enriching the interstellar medium with a- 
elements explode as core-collapse supernovae (SNe) (see, e.g. 
Ballero et al. 2007b). Many papers report the existence of the 
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Table 2. The average values of the ages, eccentricities, and Z max for the three stellar groups, along with their rats and the number of 
stars. 



Stellar groups BASTI age Padova age Eccentricity Z] 



Thin disk 




5.4±2.2 (555) 
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0.14±0.1 (549) 
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0.32±0.34 (45) 
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Fig. 5. [a/Fe] versus eccentricity and Z lnax . The symbols are the 
same as in Fig. 4. 



"knee" in [a/Fe] trends for the thick disk stars (kinematically 
selected) when [Fe/H] reaches to * -0.3 (see, e.g. Feltzing et 
al. 120031 Bensby et al. 2007). This downturn in [a/Fe], up to 
solar metallicities (Bensby et al. 2007), means that low- and 
intermediate-mass stars start to play a significant role in the 
chemical enrichment, through the explosions of SN la, which 
produce relatively little of the a-elements. This downturn is seen 
in Fig 1. As already noted, however, there is a gap between 
metal-poor and metal-rich a-enhanced stars, which casts doubt 
on the relation between the hamr and the thick disk stars. 

To study the differences and/or the similarities of the 
Galactic orbital parameters and ages between two a-enhanced 
groups, we cross-matched our sample with the Geneva- 
Copenhagen Survey (GCS) sample (Casagrande et al. 1201 II) . 
which provides the eccentricities of the orbits, a maximum ver- 
tical distance (Z max ) a star can reach above the Galactic plane, 
and the ages of about 700 of our stars. 

The mean ages of the three groups and their standard devi- 
ations are presented in Table 2. As can be seen from the table, 
the thick disk and hamr family stars have almost the same age, 
being on average older than thin disk stars by about 3 Gyr. 

Our study of the orbital parameters has uncovered a differ- 
ent connection between the three groups. Most of the hamr stars, 
such as the thin disk stars, have nearly circular orbits, close to the 
Galactic plane. Unlike them, thick disk stars move in more ec- 
centric orbits and have on average vertical distances three times 
larger than other two groups. The [a/Fe] distributions as a func- 
tion of eccentricities and Z max are shown in Fig. 5. The mean val- 
ues of the eccentricities and Z maj! for the above separated three 
stellar families are presented in Table 2. 

We can summarize the results obtained in this letter in the 
following way: our sample of stars, which consists of 1 1 12 FGK 
long-lived dwarfs, can be clearly separated into two groups ac- 
cording to the content of a elements. In turn, high-a stars can 
also be divided into two families based on the gap between 
more a-enhanced (thick disk) and less a-enhanced (hamr) stars. 



Studying these three stellar families, we see that hamr stars have 
orbits similar to the thin disk stars, but that they are similar to 
thick disk stars in terms of age. This unusual results for metal- 
rich stars assigned to both the thin and thick disks (by their kine- 
matics) was found previously in Soubiran et al. (12005b . 

In the past couple of years, many studies have focused on 
the abundance determination of the bulge and inner disk stars 
(Fulb right et al. [20051 HSU M elendez et a l.l2008l Alves-B rito 
et al. l20T0l Epstein et al. l20T0l Bensby et al. l2010al B l20TTal lbl. 
The abundance studies for the bulge dwarfs at sub-solar metal- 
licities show excellent agreement with the abundance patterns 
in the local thick disk (Melendez et al. 120081 Alves-Brito et al. 
120101 Bensby et al. 201Cb). This suggests that these two popu- 
lations have comparable chemical histories. It is also intriguing 
that at super-solar metallicities bulge stars seems to be more a- 
enhanced than local thin disk stars (Fulbright et al. 2005, 2007 
see Fig. 10 in Bensby et al. 201Cb and Fig. 10 in Alves-Brito et 
al. 120101 ). Taking into account that the high-a metal-rich stars of 
our sample are on average as old as bulge stars, and adding to 
this the aforementioned similarities between the different popu- 
lations, we can see a link between the hamr stars and stars in the 
inner disk. Stellar radial migration could give an explanation of 
this link (e.g. Haywood l200^bl Schonrich & Binney |2009al) . It 
has been proposed that metal-rich stars found in the solar vicin- 
ity may have been formed in the inner Galactic disk regions (e.g. 
Gren on [19991 Ecuvillon et al. 120071 Fam aey et al. 120571 Santos 
et al. 120081 Schonrich & Binney |2009bl ). 

Nevertheless, the origin and nature of these stars remains un- 
clear and needs to be clarified. Although the present observations 
suggest that hamr stars (high-alpha, metal rich) may have orig- 
inated from the inner disk (e.g. inner thick-disk members), they 
do not allow us to exclude the possibility that they represent a 
whole new Galactic population. More observations are needed 
to resolve this uncertainty. 

Acknowledgements. V.Zh.A. and S.G.S. acknowledge the support from the 
Funda£o para a Cilncia e Tecnologia (FCT) in the form of a grants 
SFRH/BPD/70574/2010 and SFRH/BPD/4761 1/2008, respectively. N.C.S. 
acknowledges the support of the European Research Council/European 
Community under the FP7 through a Starting Grant, as well as the support from 
(FCT), Portugal, through program Cilncia 2007. We also acknowledge support 
from FCT in the form of grant reference PTDC/CTE-AST/098528/2008. 

References 

Adibekyan, V. Zh., et al. 201 1, in preparation 

Anders, E., & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 197 
Alves-Brito, A., Melendez, J., Asplund, M., et al. 2010, A&A, 513, A35 
Bensby, T., Feltzing, S., & Lundstrom, I. 2003, A&A, 410, 527 
Bensby, T, Feltzing, S., Lundstrom, I., & Ilyin, I. 2005, A&A, 433, 185 
Bensby, T, Zenn, A. R., Oey, M. S., & Feltzing, S. 2007, ApJ, 663, L13 
Bensby, T., Alves-Brito, A., Oey, M. S., et al. 2010a, A&A, 516, L13 
Bensby, T, Feltzing, S., Johnson J.A., et al. 2010b, A&A, 512, A41 
Bensby, T, Alves-Brito, A., Oey, M. S„ e t al. 2011, ApJ, 735, L46 
Bensby, T, et al. 2011,[arXiv:1107.5606l 



4 



V. Zh. Adibekyan et al.: A new a-enhanced super-solar metallicity population 



Casagrande, L., Schonrich, R., Asplund, M., et al. 2011, A&A, 530, A138 
Dehnen, W., & Binney, J. J. 1998, MNRAS, 298, 387 
Epstein, C. R., Johnson, J. A., Dong, S., et al. 2010, ApJ, 709, 447 
Ecuvillon, A., Israelian, G., Pont, E, Santos, N. C, & Mayor, M. 2007, A&A, 
461, 171 

Famaey, B., Pont, E, Luri, X., et al. 2007, A&A, 461, 957 
Feltzing, S., Bensby, T., & Lundstrm, I. 2003, A&A, 397, LI 
Feltzing, S., & Gustafsson, B. 1998, A&AS, 129, 237 
Fuhrmann, K. 1998, A&A, 338, 161 

Fulbright, J. P., Rich, R. M, & McWilliam, A. 2005, Nucl. Phys. A, 758, 197 
Fulbright, J. P., McWilliam, A., & Rich, R. M. 2007, ApJ, 661, 1152 
Gazzano, J., de Laverny, P., Deleuil, M., et al. 2010, A&A, 523, A91 
Gazzano, J., PhD thesis, 201 1, Marseille 
Gilmore, G. & Reid, N. 1983, MNRAS, 202, 1025 
Grenon, M. 1999, Ap&SS, 265, 331 
Haywood, M. 2008a, A&A, 482, 673 
Haywood, M. 2008b, MNRAS, 388, 1175 

Kurucz, R. 1993, ATLAS9 Stellar Atmosphere Programs and 2 km/s grid. 

Kurucz CD-ROM No. 13. Cambridge, Mass.: Smithsonian Astrophysical 

Observatory, 1993., 13 
Lo Curto, G, Mayor, M., Benz, W., et al. 2010, A&A, 512, A48 
Mayor, M., Pepe, E, Queloz, D., et al. 2003, The Messenger, 114, 20 
Mayor, M., Marmier, M., Lovis, C, et al. 201 1, ArXiv e-prints 
Melendez, J., Asplund, M., Alves-Brito, A., et al. 2008, A&A, 484, L21 
Navarro, J. E, Abadi, M. G, Venn, K. A., et al. 201 1, MNRAS, 412, 1203 
Neves, V., Santos, N. C, Sousa, S. G, et al. 2009, A&A, 497, 563 
Nissen, P. E., & Schuster, W. J. 2010, A&A, 511, L10 
Petigura, E. A. & Marcy, G. W. 201 1, ApJ, 735, 41 
Ramirez, I., Allende Prieto, C, & Lambert, D. L. 2007, A&A, 465, 271 
Reddy, B. E., Lambert, D. L., & Allende Prieto, C. 2006, MNRAS, 367, 1329 
Robin A. C, Reyle C, Derriere S., Picaud S., 2003, A&A, 409, 523 
Santos, N. C, Melo, C, James, D.J., et al. 2008, A&A, 480, 889 
Santos, N. C, Mayor, M., Bonfils, X., et al. 2011, A&A, 526, 112 
Schwarzschild K., 1907, Gottingen Nachr, 614 
Schonrich, R., & Binney, J., 2009a, MNRAS, 396, 203 
Schonrich, R., & Binney, J., 2009b, MNRAS, 399, 1 145 
Soubiran, C, & Girard, P., 2005, A&A, 438, 139 
Sousa, S. G, Santos, N. C, Israelian, G, et al. 2007, A&A, 469, 783 
Sousa, S. G, Santos, N. C, Mayor, M., et al. 2008, A&A, 487, 37 
Sousa, S. G, Santos, N. C, Israelian, G, et al. 2011a, A&A, 526, 99 
Sousa, S. G, Santos, N. C, Israelian, G, et al. 2011b, A&A, 533, 141 
Sneden, C. 1973, Ph.D. Thesis, Univ. of Texas 
van Leeuwen, F. 2007, A&A, 474, 653 



5 



